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Design Device of ControLler^ 




Technical Field 



5 



The present invention relates to a design device which 



designs a multivariable controller in accordance with H 
infinity (H oo) logic. 
Background Art 

In recent years, in the feedback control field, the H 

10 oo logic has often been utilized which allows the design of 
a controller in view of an error between an actual control 
object and a numerical model of the control object. In the 
conventional control logic, when designing a control system, 
a control object model that is represented by the transfer 

15 function and the state equation is prepared, and the control 
system is designed to stabilize the model. At this time, when 
there is an insignificant error between the actual control 
object and the model, a controller designed to stabilize the 
model can also stabilize the actual control object. However, 

20 when there is a significant error between the actual control 
object and the model for some reason, the controller may not 
stabilize the actual control object. 

In the H oo logic, even if there is an error between the 
actual control object and the numerical model for use with 

2 5 the design, when the information as to the error can be 
obtained, the controller for stabilizing the actual control 
object can be designed in view of the error. It is said that 



- 1 - 



the H oo logic is more likely to give the control specification 
intuitively in designing the control system, as compared with 
the conventional control logic. For example, in a case of 
designing a control system using the conventional control 
5 logic, its design specification involved a pole in the closed 
loop system or a weight matrix of evaluation function. However, 
the physical meanings of these values were unclear, and it 
O required a lot of trial and error to make the settings. 

CO On the contrary, in the H oo logic, the control 

H 10 specif ication can be defined in accordance with the frequency 
00 responses of the closed loop system consisting of the control 

s ob j ect and the controller . The H oo logic has such an advantage, 

yj but is theoretically difficult, and has not been put into 

hj 

ifl practical use in the current situation for the reasons of 

12 15 requiring the considerable knowledge to construct the actual 
control system, and it being difficult to give the control 
specification to the objects less treatable in accordance 
with the frequency responses in the process control and so 
forth. 

20 Disclosure of the Invention 

A control obj ect model used in designing a multivariable 
control system has a variety of magnitudes of error components 
ranging from the manipulated variables to controlled 
variables. In this way, there are various error components 

25 of the model to the manipulated variables. In the H oo logic, 
since the controller is designed with reference to the 
component with a large gain of error, the component with a 
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small gain of error is prone to have a very conservative 
response, or to be excessively stable. In some cases, it may 
be required to adjust the weight of control for each 
controlled variable to avoid the interference between 
controlled variables. Thus, it is proposed to introduce a 
manipulated variable weight called a scaling matrix T in order 
to make the error magnitudes of the control object models 
consistent and weight the control for the controlled 
variables. However, in the conventional design method, no 
general solution was not established to determine the scaling 
matrix T, resulting in the problem that the scaling matrix 
T was difficult to suitably choose. Since it was difficult 
to determine the scaling matrix T, there was another problem 
that the H oo logic was difficult to utilize in the design 
of the multivariable controller. 

Also, in the H oo logic, it is required to determine a 
frequency weight called a sensitivity weight W s to determine 
the set value followup characteristic of the closed loop 
system. However, since the H oo logic is a design method in 
the frequency domain, there was the problem that the design 
is easy in the control of the mechanical system, but is hard 
in the control system untreatable in the frequency domain 
such as the process control, whereby it is difficult to 
suitably select the sensitivity weight W s . Also, since it is 
difficult to give the control specification in the frequency 
domain, and hard to determine the sensitivity weight W s , there 
was the problem that the H oo logic is difficult to utilize 



in the design of the controller used in the process control 
field. 

As described above, there was conventionally the problem 
that the H oo logic was hard to utilize in the design of the 
controller . 

The present invention has been made in order to solve 
these problems, and has as its object to provide a design 
device that can easily design a controller in accordance with 
the H co logic. 

A design device of a controller according to the present 
invention comprises storage means for storing generalized 
plants, parameter calculating means for calculating the 
parameters of the components of the generalized plants in 
accordance with the response characteristic of a control 
object model or the response characteristic of a closed loop 
system consisting of the control object model and the 
controller, and controller calculation means for calculating 
the parameters of the controller by applying the parameters 
to the generalized plants stored in the storage means. 

In the design device of the controller according to one 
arrangement of the invention, the generalized plants have 
the control object models, and manipulated variable weight 
adjusting means for adjusting the input of manipulated 
variables into the control object models, which is provided 
in the former stage of the control object models, the 
parameter calculating means comprises frequency response 
calculation means for calculating the frequency responses 



of the control object models, and scaling matrix calculation 
means for calculating a scaling matrix T for determining the 
weighting of the manipulated variables with the manipulated 
variable weight adjusting means in accordance with the 
5 frequency responses of the control object models so that the 
respective gains of the control object models are consistent, 
and the controller calculation means calculates the 
parameters of the controller by applying the scaling matrix 
T to the manipulated variable weight adjusting means of the 
10 generalized plants stored in the storage means. 

Also, in the design device of the controller according 
to one arrangement of the invention, the generalized plants 
have a first control object model for the manipulated 
variables, a second control object model for the disturbance, 
15 and manipulated variable weight adjusting means for 

adjusting the input of manipulated variables into the first 
control object model, which is provided in the former stage 
of the first control object model, the parameter calculating 
means comprises frequency response calculation means for 
20 calculating the frequency responses of the first control 
object model and the second control object model, and scaling 
matrix calculation means for calculating a scaling matrix 
T for determining the weighting of the manipulated variables 
with the manipulated variable weight adjusting means in 
25 accordance with the frequency responses of the first and 
second control object models so that the respective gains 
of the first control object model are consistent with the 
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maximum value of the gains of the second control obj ect model, 
and the controller calculation means calculates the 
parameters of the controller by applying the scaling matrix 
T to the manipulated variable weight adjusting means of the 
5 generalized plants stored in the storage means. 

Also, in the design device of the controller according 
to one arrangement of the invention, the generalized plants 
stored in the storage means have control variable weight 
adjusting means for adjusting the controlled variables 
10 inside a closed loop system consisting of the manipulated 
variable weight adjusting means, the control object model 
for the manipulated variables and the controller, and the 
design device has setting means for setting a weight matrix 
S for determining the weighting of the controlled variables 
15 with the control variable weight adjusting means. 

Also, in the design device of the controller according 
to one arrangement of the invention, the generalized plants 
stored in the storage means have control variable weight 
adjusting means for adjusting the controlled variables in 
2 0 the former or latter stage of frequency sensitivity weight 
adjusting means to determine the set value followup 
characteristic of a closed loop system consisting of the 
manipulated variable weight adjusting means, the control 
object models for the manipulated variables and the 
25 controller, and the design device has setting means for 

setting a weight matrix S for determining the weighting of 
the controlled variables with the control variable weight 
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adjusting means. 

In the design device of the controller according to one 
arrangement of the invention, the parameter calculating 
means has setting means for setting the transient response 
5 characteristic of the closed loop system and frequency 
sensitivity weight calculation means for calculating the 
frequency sensitivity weight to determine the set value 
^ followup characteristic of the closed loop system in 

ft accordance with the transient response characteristic of the 

"2 10 closed loop system, and the controller calculation means 

K calculates the parameters of the controller by applying the 

PJ 

frequency sensitivity weight to the generalized plants 
j3 stored in the storage means. 

m In the design device of the controller according to one 

M= 15 arrangement of the invention, the frequency sensitivity 
weight calculation means calculates the frequency 
sensitivity weight in accordance with the transient response 
characteristic of the closed loop system, and a design index 
that the H 00 norm of a transfer function of the closed loop 
20 system from the set value to the deviation, multiplied by 
the frequency sensitivity weight, is less than 1. 

In the design device of the controller according to one 
arrangement of the invention, the setting means approximates 
the transient response characteristic of the closed loop 
25 system with a first-order lag characteristic. 

In the design device of the controller according to one 
arrangement of the invention, the setting means approximates 
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the transient response characteristic of the closed loop 
system with a second-order system characteristic. 
Brief Description of the Drawings 

Fig. 1 is a block diagram showing the configuration of 
5 a design device of a controller in a first embodiment of the 
present invention; 

Fig. 2 is a block diagram showing the configuration of 
a model of an actual control object as represented in the 
formulae; 

10 Fig. 3 is a block diagram showing the configuration of 

the conventional generalized plants; 

Fig. 4 is a block diagram showing the configuration of 
the generalized plants for use with the design device of the 
invention; 

15 Fig. 5 is a block diagram showing the configuration of 

a robust control system with a controller added to the 

generalized plants of Fig. 4; 

Fig. 6 is a block diagram showing the configuration of 

an actual controller containing the controller designed 
20 using the design device in the first embodiment of the 

invention; 

Fig. 7 is a block diagram showing an additive error in 
a numerical model of the control object; 

Fig. 8 is a block diagram showing the configuration of 
25 a closed loop system used in determining the sensitivity 
weight; 

Fig. 9 is a graph for explaining the operation of a 
- 8 - 



scaling matrix in the first embodiment of the invention; 

Fig. 10 is a graph for explaining the operation of a 
scaling matrix in a second embodiment of the invention; 

Fig. 11 is a block diagram showing the configuration 
5 of a design device of a controller in a third embodiment of 
the present invention; 

Fig. 12 is a block diagram showing the configuration 
2 of the generalized plants in the third embodiment of the 

%j\ present invention; 

JJ io Fig. 13 is a block diagram showing the configuration 

SI of an actual controller containing the controller designed 

P using the design device in the third embodiment of the 

n\ invention; 

~LJ Fig. 14 is a block diagram showing the configuration 

15 of the generalized plants in a fourth embodiment of the 
present invention; 

Fig. Ij3_is a block diagram showing the configuration 
of a design device of a controller in a fifth embodiment of 
the present invention; 
20 Fig. 16 is a block diagram showing the configuration 

of a model of an actual control object as represented in the 
formulae; 

Fig. 17 is a block diagram showing the configuration 
of the conventional generalized plants; 
25 Fig. 18 is a block diagram showing the configuration 

of the generalized plants for use with the design device of 
the invention; 
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Fig. 19 is a block diagram showing the configuration 
of a robust control system with a controller added to the 
generalized plants of Fig. 18; 

Fig. 20 is a block diagram showing the configuration 
of an actual~"controller containing the controller designed 
using the design device of the invention; 

Fig. 21 is a block showing an additive error in a 
numerical model of the control object; 

Fig. 22 is a block diagram showing the configuration 
of a closed loop system used in determining the sensitivity 
weight; and 

Fig. 23 is a graph showing the time response 
characteristic when the closed loop system is approximated 
by a second-order system. 
Best Mode of Carrying Out the Invention 
First embodiment 

The present invention will be described below in detail 
by way of embodiment with reference to the accompanying 
drawings. Fig. 1 is a block diagram showing the configuration 
of a design device of a controller in a first embodiment of 
the invention. The design device of Fig. 1 comprises a control 
object model input unit 1 for inputting the parameters of 
a control object model, a control object model registration 
unit 2 for registering the model parameters in a memory unit 
described later, the memory unit 3 for storing formulae of 
generalized plants and formulae of control object models 
which are parts of the generalized plants, a frequency 



response calculation unit 4 for calculating the frequency 
responses of the control object models, a scaling matrix 
calculation unit 5 for calculating a scaling matrix T for 
making consistent the magnitudes of errors of the control 

5 object models so that the respective gains of the control 
object models are consistent with the maximum value of the 
gains of the control object models, and a controller 
calculation unit 6 for calculating the parameters of a 
controller by applying the scaling matrix T to the generalized 

10 plants stored in the memory unit 3. 

An algorithm for designing the controller in accordance 
with the H oo logic is effected based on the generalized plants 
as represented using the control objects. Therefore, the 
generalized plants will be first described. Fig. 2 is a block 

15 diagram showing the configuration of a model of an actual 
control object as represented in the formulae. A numerical 
model of the control object as shown in Fig. 2 is composed 
of a first control object model 11 for the manipulated 
variable u and a second control object model 12 for the 

20 disturbance w. Reference sign Pu denotes a transfer function 
of the model 11 and reference sign Pw denotes a transfer 
function of the model 12. The models 11, 12 are obtained as 
a result of the model identification using the data obtained 
by a step response test for the actual control object. The 

25 controlled variable y that is the output of control object 
is a sum of the outputs from the models 11 and 12. 

Fig. 3 shows the configuration of the conventional 
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generalized plants containing such numerical model of 
control object. The generalized plants involve providing a 
frequency weight called a sensitivity weight W s used to 
determine the set value followup characteristic and a 
frequency weight called a complementary sensitivity weight 
W t used to determine the robust stability, and introducing 
the set value r, and the outputs zl, z2 in addition to the 
input (manipulated variable) u of control object, input 
(disturbance) w, and output (controlled variable) y, to 
realize both the set value followup characteristic and the 
robust stability as shown in Fig. 3. Deviation e (=y-r) is 
an observation quantity or the input into the controller (not 
shown) . Reference numeral 13 denotes a block (frequency 
sensitivity weight adjusting means) representing the 
sensitivity weight W s , and reference sign Zi denotes an output 
for effecting the evaluation of the set value followup 
characteristic. Also, reference numeral 14 denotes a block 
representing the complementary sensitivity weight W t , and 
reference sign Z 2 denotes an output for effecting the 
evaluation of the robust stability. 

Conventionally, in the generalized plants as shown in 
Fig. 3, the complementary sensitivity weight W t is determined 
by estimating the uncertainty of the models on the basis of 
the numerical models of control objects, and the sensitivity 
weight W s is determined by directly specifying the frequency 
characteristic in view of the f ollowupability to the set value 
r, whereby the parameters of the controller are determined 



through the y iteration. However, since employing the 
generalized plant of Fig. 3, the controller is designed on 
the basis of a larger gain due to differences in the gains 
for the plant outputs of the manipulated variables, the 
obtained controller is likely to be very conservative, or 
excessively stable. Also, since the set value followup 
characteristic and the disturbance response characteristic 
are usually reciprocal, it is preferred to design with the 
weight according to the purpose rather than with the same 
weight. Further, in the generalized plants of Fig. 3 
containing no integral element, some steady-state deviation 
arises. If the sensitivity weight W 3 is provided with the 
integral characteristic, the controller can have the 
integral characteristic, but the generalized plants become 
unstable, not leading to the normal H oo problem. 

Therefore, in this embodiment, the generalized plants 
as shown in Fig. 4 are considered. In Fig. 4, reference sign 
M denotes a scaling matrix for adjusting the influence of 
disturbance w on the controlled variable y, reference sign 
T denotes a scaling matrix for making the magnitude of error 
of the control object model consistent, and reference sign 
a -1 1 denotes a weight for providing the controller with the 
integral characteristic to eliminate the steady-state 
deviation. Herein, <x(s) is defined as a(s)= s/ (s+a) . Where 
s is a Laplace operator, and a (>0) is any real number. 
Reference numeral 15 denotes a block representing the scaling 
matrix M, reference numeral 16 denotes a block (manipulated 
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variable weight adjusting means) representing the scaling 
matrix T, and reference numeral 17 denotes a block 
representing the weight a -1 1. Deviation e 2 is deviation e 
multiplied by the weight a _1 I, and input into the controller . 
Fig. 5 shows the configuration of a robust control system 
in which a controller K is added to the above generalized 
plant. In Fig. 5, reference numeral 18 denotes a block 
representing the controller K. 

The design device of the controller in this embodiment 
is aimed at determining the parameters of the controller K 
such that the controlled variable y that is the output of 
control object follows the set value r, the influence of 
disturbance w is removed, and the control obj ect is stabilized 
even if it is fluctuated or there is an error in the model 
of control object. The H oo problem can be regarded as the 
problem of reducing the H oo norm (gain) of the transfer 
function from (r, w) to (z lr z 2 ) . That is, the set value 
followup characteristic, the robust stability and the 
disturbance suppression may be considered in the following 
way. 

(A) Set value followup characteristic: If the H co norm 
(gain) of the transfer function from the set value r to 
deviation e (more correctly, transfer function with the set 
value r multiplied by the frequency weight a -1 W s from r to 
zi) is reduced, the deviation e can be decreased, so that the 
set value followup characteristic can be made better . Herein, 
a _1 W s is the frequency weight for restricting the followup 



band (e.g., followup only in the low band) . 

(B) Robust stability: There is an error between the 
actual control object and its model due to characteristic 
variations of the control object or the error at the time 
of modeling. The maximum value of the error from the 
identified model is estimated as A(s), and the controller 
K is designed such that the H oo norm from the set value r 
to z 2 is smaller than or equal to 1, employing the 
complementary sensitivity weight W t (s) such- as 
| A (jo) I < I W t ( j co ) | for this error, whereby the robust 
stabilization can be achieved. 

(C) Disturbance suppression: If the H oo norm (gain) of 
the transfer function from disturbance w to deviation e (more 
correctly, transfer function with the disturbance w 
multiplied by the frequency weight a -1 W s from w to Zi) is 
reduced, the deviation e can be decreased even if the 
disturbance w enters, whereby the disturbance suppression 
can be ameliorated. 

Next, it is supposed that the state space representation 
of the generalized plant as shown in Fig. 4 is given by: 

• • • (1) 

y = C p x p + D pl Mw+D p2 Tu • • • (2) 

In the above expressions (1) and (2) , x p is the quantity 
25 of state, and A p , B pl , B p2 , C p , D pl and D p2 are the parameters 
of the numerical models 11, 12 of control objects. From the 
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expression (2), the deviation e can be obtained in the 
following expression. 



= y- r = 



, +D n1 Mw + D n7 Tu- r * * * < 3 > 



With the configuration of the generalized plant as shown 
in Fig. 4, the outputs z/ and z 2 ' can be defined in the 
following expressions. 

z,' =e 2 • • • (4) 

The frequency weight for providing the controller Kwith 
the integral characteristic can be defined in the following 
expression, using the expression (3) . 

X a =A a x a + B a e 

= A A + B a C p x p + B B D pl Mw + B a D p2 Tu-B a r • • • (6) 



" 5 



0 



e2 = C a x a +D a e 

= C a x a + D a C p x p + D a D pl M w + D a D p2 T u - D a r 



(7) 



In the expressions (6) and (7), x« is the quantity of 
state of a -1 1, and A«, B«, C a and D„ are the parameters of a _1 I . 
Arranging the above expressions and representing them in a 
state space, the following three expressions can be obtained. 



BpiM B p2 T 



A[ x p1 = [ a p °1PS 

dt[x a J [B a C p A a j[x a J r [-B a B a D pl M B a D p2 T 

o cpC oiH 

= |P«Cp C a ^p| + [-D a D a D pl M D a D p2 T|w 



-D a D a D pl M D a D p2 Tl 
0 0 I 



(8) 



(9) 
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Representing the expressions (8) , (9) and (10) in 
Doile' s notation, the following expression can be obtained. 



G(s)= 



A Bi B 2 



C 2 D 2 



Parameter A can be represented as follows: 



A=f A P ° 
B„C p A a 



(1 1) 



(1 2) 



Parameters Bi and B 2 can be represented as follows: 



0 B pl M 



Bo 



B n2 T 



(13) 



Also, parameters Ci and C 2 can be represented as follows: 
TD„C n C„ 



= [°o Cp o"] ■ c ^ c » C »J 



(14) 



And the parameters Dn, D 12 , D 21 , D 22 can be represented 
as follows: 



D 2 i=L-Da D a D pl Mj , D^DJ^T ... (15) 

The sensitivity weight W s and the complementary 
sensitivity weight W t are designed, and multiplied by the 
output parts of the expression (11) . Through y iteration, 
the controller K is obtained in the state space representation . 
Herein, the output parts of the expression (11) signify the 
parts corresponding to the outputs Zi' , z 2 ' in Fig. 4. Hence, 
an output equation of the parameters Ci, Dn, D i2 of the 
expression (11) may be multiplied by a diagonal matrix Q as 
represented in the following expression having the diagonal 
elements of sensitivity weight W s and complementary 
sensitivity weight W t from the left side. Thus, the parameters 
of the controller K can be calculated. 



[0 W t J 
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The controller K is a solution of the H <x> control problem 
with the generalized plants, and the actual controller 
mounted on the plants such as a distillation tower is the 
controller K multiplied by weight a _1 I and scaling matrix T, 
as shown in Fig. 6. 

Next, a method for determining the complementary 
sensitivity weight W t in this embodiment will be described 
below. The control object is varied in the characteristics, 
depending on the driving conditions. Normally, the control 
design is made based on a certain model, but in the robust 



control design, the variation of control object and the 
magnitude of error of modeling are contained beforehand in 
the control design, and the controller is designed to be 
stable without much deterioration in the control performance 
even if there is any variation or error. Fig. 7 shows an 
additive error for the model 11 of control object. In Fig. 
7, reference numeral 19 denotes a block representing the 
additive error A. In the robust control design, a variation 
in the characteristic of the control object owing to the 
driving conditions and a model error due to lower dimension 
of the model 11 are represented as the additive error A as 
shown in Fig. 7. If the characteristics of the control object 
are deviated from the model 11 due to this additive error 
A, the controller is designed such that the controller output 
may be stable. To this end, the complementary sensitivity 
weight W t may be determined to cover the additive error A. 
The general expression of this complementary sensitivity 
weight W t is shown in the following expression. Since the 
change of the model 12 is not related with the stability of 
the system, it is supposed that the model 11 alone is varied 
in designing the controller. 
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In this embodiment, for the additive error A with scaled 
magnitude, employing the scaling matrix T, the element of 



the complementary sensitivity weight W t is the maximum value 
Gmax of the gains of error A multiplied by a safety factor 
5 (5 is equal to 1, for example) . Namely, the elements 
(weights) W ti , W t2 , W t3 / •••/ w tN are defined as follows. 
-W ll = Wt2 = W t3 = WtN= (1 + 5) Gmax • • • (18) 

This embodiment is involved with the multivariable 
control system, and assuming that the number of manipulated 
variables u is N (N is a positive integer) , the complementary 
sensitivity weight W t is NxN matrix. Where W tN is the weight 
for the N-th manipulated variable u N . 

Next, a method for determining the sensitivity weight 
W s in this embodiment will be described below. First of all, 
a general expression of the sensitivity weight W s is shown 
in the following expression. 
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Assuming that the number of controlled variables y is 
L (L is a positive integer) , the sensitivity weight W s is LxL 
matrix. Element W sL of the sensitivity weight W s is the weight 
for the L-th controlled variable y L . In order to determine 
the sensitivity weight W s , a closed loop system of Fig. 8 in 
which the robust control system of Fig. 5 is simplified is 
considered. In Fig. 8, reference numeral 11a denotes a block 
representing a numerical model P of control object, and 
reference numeral 13a denotes a block representing the 
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frequency weight W a ' . Supposing that the sensitivity function 
indicating the control performance mainly regarding the 
quick-response property such as the set value 
followupability or disturbance suppression is S(s), the 
5 smaller gain |S(j©)| of the sensitivity function S (s) is 
preferable because the model variations have less effect on 
the set value response. If the control specification at each 
frequency is given by S spe c(G>)/ the following condition 
concerning the sensitivity function S(s) can be obtained. 
10 |S(jo>)|<S spec (a>);Va> • • • (2 0) 

Where V© means that the expression (2 0) holds for all 
frequencies ©. Employing this sensitivity function S (s) , the 
design index of the controller K in view of the set value 
followup characteristic is as follows: 

15 |w sL '(s)S(s)||„<l • ■ • (2D 

The frequency weight W sL ' (s) is W sL (s) multiplied by 
cc -1 (s), and is defined as follows: 

W sL '(s)= a - 1 (s)W sL (s) • • • (2 2) 

The expression (21) indicates that the H oo norm of the 
20 transfer function of the closed loop system of Fig. 8 from 
the set value r to the deviation e (more correctly, the 
transfer function with the set value r multiplied by the 
frequency weight of 1 (s) W sL (s) from r to z a ) is less than 1. 
By setting the weight W sL (s) to satisfy this expression (21) , 
25 the controller K can be designed in view of the set value 
followup characteristic. 
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Next, a method for determining the scaling matrix M will 
be described below. A general expression of the scaling matrix 
M is shown in the following expression. 
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Assuming that the number of disturbance w is J (J is 
a positive integer), the scaling matrix M is JxJ matrix. 
Element Mjof the scaling matrix M is the weight for the J-th 
disturbance w a , with the initial value being 1. Each element 
Mj is an adjustment parameter for determining the disturbance 
suppression performance by adjusting the influence of each 
disturbance Wj on the controlled variable y. That is, when 
the suppression of the particular disturbance w is desired 
to be intensified, the element Mj concerning this disturbance 
w is made larger than 1. 

Next, a method for determining the scaling matrix T in 
this embodiment will be described below. A general expression 
of the scaling matrix T is shown in the following expression. 



T = 



Ti 0 0 
0 T 2 0 
0 0 T 3 



(24) 



0 0 0 - T N j 

Assuming that the number of manipulated variables u is 
N (N is a positive integer) , the scaling matrix T is NxN matrix . 
Element T N of the scaling matrix T is the weight for the N-th 
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manipulated variable u N . Each element T N is determined such 
that the magnitudes of the gains of the control object model 
11 are as equal as possible. More specifically, each element 
T N is determined as follows: 




• - • (2 5) 



In the expression (25) , G yLuN is the transfer function 
of the control object model 11 of Fig. 4 from the N-th 
manipulated variable u N to the L-th controlled variable y L , 
and llGyLuwII 00 is the H oo norm (gain) of the same transfer function, 
max (||G yLul ||°o, ||G yLu2 |K ||G yLuN ||°o) means to select the maximum 
value from among the H qo norms ||G yLu i||oo, ||G yLu2 ||°o, and ||G yLuN |la>. 
To obtain the H oo norm ||G y lun||°°, the model 11 as presented 
in the state equation representation may be transformed into 
the transfer function representation given by the following 
expression, and the gain of each frequency may be calculated 
from this transfer function. 
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Thereby, the H oo norm ||G yl uN||«> can be obtained for each 
of the manipulated variables u and the controlled variables 
y, and the element T N of the scaling matrix T can be obtained 
from the expression (25) . 
5 Referring now to Fig. 9, the operation of the scaling 

matrix T will be described below. Fig. 9A shows the gain 
characteristic of the control object model 11 (frequency 
response characteristic of the model 11) . In Fig. 9, three 
sorts of gain characteristic are only shown to simplify the 
10 description, but if the number of manipulated variables u 
is N and the number of controlled variables y is L, NxL sorts 
of gains exist. As shown in Fig. 9A, in the case where there 
is no scaling matrix T, it will be found that the gains of 
the control object model 11 are inconsistent. In general, 
15 if the gains of the control object model are inconsistent, 
the magnitudes of the errors of the control object model are 
correspondingly inconsistent. Since the complementary 
sensitivity weight W t is determined to cover the additive 
error A, as previously described, it follows that the 
20 controller is designed on the basis of the model with larger 
error, whereby the controller obtained is likely to be very 
conservative, or excessively stable. 

Thus, the magnitudes of the gains are made consistent, 
employing the scaling matrix T. Fig. 9B shows the gain 
25 characteristic of the control object model 11 in the case 
where the scaling matrix T of this embodiment is provided. 
llGyumaxll 00 is the maximum value of the gains of the model 11. 
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As will be clear from Fig. 9B, the method for determining 
the scaling matrix T of this embodiment as shown in the 
expressions (24) and (25) involves determining the scaling 
matrix T such that the gains are consistent with the gain 
5 maximum value ||G yuma x||°o (more correctly, the approximation of 
the gain maximum value) of the model 11. 

Referring now to Fig. 1, the above operation will be 
described below. The parameters of the control object model 
11 are set in the control object model input unit 1 by the 
10 user, of the design device. The control object model 

registration unit 2 registers the parameters input from the 
control object model input unit 1 into the formulae of the 
control object model stored beforehand in the memory unit 
3. The control object model input unit 1 and the control object 
15 model registration unit 2 constitute model setting means for 
setting the control object model. The memory unit 3 stores 
the formulae of the generalized plants of Fig. 4 and the 
formulae of the control object models that are parts of the 
generalized plants, as described in the expressions (1) to 
20 (15) . The frequency response calculation unit 4 transforms 
the model 11 as presented in the state equation representation 
which is registered in the memory unit 3 into the transfer 
function representation, and calculates the gain at each 
frequency from the transfer function. Subsequently, the 
25 scaling matrix calculation unit 5 calculates the scaling 
matrix T, employing the expressions (24) and (25), on the 
basis of the gains calculated in the frequency response 
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calculation unit 4, and outputs it to the controller 
calculation unit 6. The controller calculation unit 6 
registers the scaling matrix T in the formulae of the 
generalized plants stored in the memory unit 3, and through 
y iteration, calculates the parameters of the controller K. 
At this time, the complementary sensitivity weight W t , the 
sensitivity weight W s and the scaling matrix M are preset in 
the generalized plants in the memory unit 3. In this way, 
the controller K can be designed. 

As previously described, the conventional method has 
no general way of determining the scaling matrix T, and 
empirically determined the scaling matrix T. On the contrary, 
in this embodiment, the scaling matrix T is calculated, in 
accordance with the frequency responses of the control object 
models 11, such that the respective gains of the control 
object models are consistent with the gain maximum value (more 
correctly, the approximation of the gain maximum value) of 
the control object model 11. Hence, the scaling matrix T can 
be easily determined. In this way, the multivariable 
controller can be easily designed in accordance with the H 
oo logic which is superior in the set value followup 
characteristic, and capable of stabilizing the controller 
even when the control objects are varied or the control object 
models 11 have some error. Consequently, the multivariable 
control system can be easily designed in view of the 
variations of the control objects and the uncertainty of the 
numerical models. Also, the multivariable controller can be 



realized utilizing the features of the H oo control with 
smaller computational load during execution of the control, 
and with which the small-scale control system can be 
implemented. 

In this embodiment, the scaling matrix T is determined 
such that the respective gains of the control object models 
are consistent with the gain maximum value (more correctly, 
the approximation of the maximum value) of the control object 
model 11, but the scaling matrix T may be determined such 
that the respective gains are consistent with the gain minimum 
value or gain average value of the control object model 11. 
To make the respective gains consistent with the gain minimum 
value (more correctly, the approximation of the minimum 
value) , max in the expression (25) may be replaced with min 
for selecting the minimum value of HGyLuill 00 /- llGyLuzll 00 , -/ and 
||G y lun||°°. And to make the respective gains consistent with the 
gain average value (more correctly, the approximation of the 
average value) , max in the expression (25) may be replaced 
with E for acquiring the average value of ||G y Luill°o, || G yLU 2 II 
and HGyLuNllco. 

Second embodiment 

In the first embodiment, the disturbance w is not 
considered, but the control object model for the disturbance 
w may be obtained. Thus, in this second embodiment, a method 
for determining the scaling matrix T in view of the influence 
of disturbance w in such a case will be described below. In 
this second embodiment, the general expression of the scaling 
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matrix T can be given by the expression (24) in the same manner 
as in the first embodiment. 

And in this second embodiment, each element T N of the 
scaling matrix T is determined in the following expression. 




• • • (2 7) 



In the expression (27), G yL „j is the transfer function 
of the control object model 12 of Fig. 4 from the J-th 
disturbance Wj to the L-th controlled variable y L , and ||G yLwJ ||oo 
10 is the H oo norm (gain) of the same transfer function, max 
{||G yLwl ||oo, ||G YLw2 ||oo, ||G yLwJ ||qo) means to select the maximum 

value from the H oo norms ||G yLw i||oo, ||G yLw2 ||°o, and ||G yLwJ ||oo. 
To obtain the H oo norm ||G yL wj||°°, the model 12 as presented 
in the state equation representation maybe transformed into 
15 the transfer function representation as in the following 
expression, and the gain of each frequency may be calculated 
from this transfer function. 
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In the expression (28) , y wL is the output of the control 
object model 12 for the disturbance w. Thereby, the H <x> norm 
||G yLwJ ||oo can be obtained for each of the disturbance w and the 
5 controlled variables y, and the element T N of the scaling 
matrix T can be obtained from the expression (27) . 

Referring now to Fig. 10, the operation of the scaling 
matrix T will be described below. Fig. 10A shows the gain 
characteristic of the control object model 12 (frequency 
10 response characteristic of the model 12) . In Fig. 10A, three 
sorts of gain characteristic are only shown to simplify the 
description, but if the number of disturbance w is J and the 
number of controlled variables y is L, JxL sorts of gain exist . 
llGy^axlloo is the maximum value of the gains of the model 12. 
15 On the other hand, Fig. 10B shows the gain characteristic 

of the control object model 11. As shown in Fig. 10B, it will 
be found that the gain maximum value of the control object 
model 12 and the gains of the control object model 11 are 
inconsistent . 

20 Fig. 10C shows the gain characteristic of the control 

object model 11 in the case where the scaling matrix T of 
this embodiment is provided. As will be clear from Fig. 10C, 
the method for determining the scaling matrix T of this 
embodiment as shown in the expressions (24) and (27) involves 
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determining the scaling matrix T such that the respective 
gains of the model 11 are consistent with the gain maximum 
value HGy^axlK (more correctly, the approximation of the gain 
maximum value) of the model 12. 

The scaling matrix T for the manipulated variable u is 
included in the closed loop system in mounting the controller . 
Accordingly, there is a meaning in making the magnitudes of 
the gains of the model 11 consistent, and it is not necessarily 
of importance with which the gain from the manipulated 
variable u to the controlled variable y is consistent. The 
previous first embodiment illustrates one instance with 
which the respective gains of the control object models are 
consistent. On the contrary, in this second embodiment, 
because the input of disturbance is considered, it is 
necessary that the influence of the input disturbance w is 
suppressed by the manipulated variable u from the viewpoint 
of suppressing the disturbance w. Thus, in this second 
embodiment, to cope with the worst situation, the scaling 
matrix T is determined such that the gains of the model 11 
are consistent with the gain maximum value HGywmaxll 00 (more 
correctly, the approximation of the gain maximum value) of 
the model 12. 

In this second embodiment, the configuration of the 
design device is almost equivalent to that of the first 
embodiment. Thus, the operation of the design device of this 
second embodiment will be described below with reference to 
Fig. 1. 



The parameters of the control object model (model 11, 
12 in this second embodiment) are set in the control object 
model input unit 1 by the user of the design device. The 
control object model registration unit 2 registers the 
parameters input from the control object model input unit 
1 into the formulae of the control object model stored 
beforehand in the memory unit 3. The frequency response 
calculation unit 3 transforms the model 11, 12 as presented 
in the state equation representation which is registered in 
the memory unit 3 into the transfer function representation, 
and calculates the gain at each frequency from the transfer 
function. Subsequently, the scaling matrix calculation unit 
5 calculates the scaling matrix T, employing the expressions 
(24) and (27), on the basis of the gain calculated in the 
frequency response calculation unit 4, and outputs it to the 
controller calculation unit 6. The operation of the 
controller calculation unit 6 is exactly the same as in the 
first embodiment. In this way, the controller K can be 
designed. 

As described above, in this second embodiment, the 
scaling matrix T is calculated, on the basis of the frequency 
response of the first and second control object models 11, 
12, such that the respective gains of the first control object 
model 11 are consistent with the maximum value of the gains 
of the second control object model 12. Hence, the scaling 
matrix T can be easily determined. Thereby, the multivariable 
controller can be easily designed in accordance with the H 



oo logic which is superior in the set value followup 
characteristic and disturbance suppression, and capable of 
stabilizing even when the control object is varied or the 
control object model has some error. 
Third embodiment 

In the first and second embodiments, the magnitudes of 
the gains of the model are made consistent by the scaling 
matrix T, so that the weights of control for the controlled 
variables are even. However, in practice, there is the problem 
that the controlled variables y interfere with each other 
in the control, and the control is unstable, whereby it may 
be sometimes required to coordinate the weights of control 
for the controlled variables. Thus, in this second embodiment , 
a weight matrix S is introduced to weight the controlled 
variables y directly. 

Fig. 11 is a block diagram showing the configuration 
of a design device of a controller in a third embodiment of 
the present invention. Fig. 12 is a block diagram showing 
the configuration of the generalized plant in this third 
embodiment. The design device of Fig. 11 is the design device 
of the first or second embodiment as shown in Fig. 1, 
comprising additionally a control variable weight input unit 
7 for inputting the weight for the controlled variable y, 
a control variable weight registration unit 8 for registering 
the control variable weight in the device, and a weight matrix 
calculation unit 9 for calculating the weight matrix S based 
on the control variable weight. Also, the generalized plant 



- 32 - 



of Fig. 12 is the generalized plant of the first or second 
embodiment as shown in Fig. 4, with a block (control variable 
weight adjusting means) 20 representing the weight matrix 
S. In this embodiment, the control variable weight adjusting 
5 means 20 (weight matrix S) is provided inside a closed loop 
system composed of the manipulated variable weight adjusting 
means 16 (scaling matrix T) , the control object model 11 and 
the controller K. A general expression of the weight matrix 
S is shown in the following expression. 



pi 


0 


0 • 


• o • 


0 


s 2 


0 • 


• 0 


0 


0 


s 3 • 


■ 0 


0 


0 


0 • 


•• s L 



Assuming that the number of controlled variables y is 
L (L is a positive integer) , the weight matrix S is LxL matrix. 
Element S L of the weight matrix S is the weight for the L-th 
15 controlled variable y L . Each element S L is determined in the 
following expression. 



SL maxiW yl ,W y2> -,W yL ) ' ' ' (3 ° 

In the expression (30) , W yL is the control variable weight 
20 for the L-th controlled variable y L . In this way, each 

controlled variable can be directly weighted by the weight 
matrix S . 

Next, in this third embodiment, the configuration of 
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the generalized plant is modified as shown in Fig. 12, so 
that the following expression can hold. 



e . = S e 



(31) 



From the expressions (31) and (3), the expressions (6) 
and (7) can be rewritten in the following expression. 



*a = A a x a + B a e l 

= A a x a + B a S C p x p + B a S D pl Mw + B a S D p2 T u - B a S r 



(3 2) 



e2 = C a x a +D a e 

^a x a ~^-D a S ^p A p 



C„x„ + D„S C n x D + D a S D pl Mw + D a S D p2 T u - D a S r 



• • • (3 3) 

Thereby, the expressions (8), (9) and (10) can be 
rewritten in the following expression. 

dfxpl [ A p Olrx p | J 0 B pl M B 2 T 1 

dT[x a J [B a SC p A a J[x a J [-B a S B a SD pl M B a SD p2 T] 

• • • (34) 



D«SC p C a -irx p 1 f-D a S D a SD pl M D a S D p2 Tl 



(3 5) 



e 2 = [D a SC p C a ]|^j + [-D a S D a SD pl M D a S D p2 t|w| 

• • • (3 6) 

Representing the expressions (34), (35) and (36) in the 
Doile's notation as in the expression (11), the parameter 



A in the expression (11) can be represented as follows: 



The parameters Bi and B 2 can be represented in the 
following expression. 



B- = f B P 2T 1 • (38) 

j pl Mj ' ^ [B a SD p2 Tj 

The parameters Ci and C 2 can be represented in the 
following expression. 

Ci = ^D a SC p C a j ^ C 2 =|D a SC p C a J ... (39) 

The parameters D U/ D 12 , D 2 i and D 22 can be represented 
in the following expression. 



Dn = J-D a S D a SD pl Mj Di2 = |D„SD p2 Tj 

D 2 l=L- D aS D a SD pl Mj . D 22 =D a SD p2 T ... (40) 

The sensitivity weight W s and the complementary 
sensitivity weight W t are designed in the same manner as in 
the first embodiment, and multiplied by the output part of 
the expression (11) , and through y iteration, the controller 
K is obtained in the state space representation. The actual 
controller mounted in the plants such as distillation tower 
is the controller K multiplied by the weight matrix S, the 
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weight a -1 1 and the scaling matrix T, as shown in Fig. 13. 

Referring now to Fig. 11, the above operation will be 
described below. The operation of the control object model 
input unit 1, the control object model registration unit 2, 
the frequency response calculation unit 4 and the scaling 
matrix calculation unit 5 is exactly the same as in the first 
or second embodiment. A memory unit 3a stores the formulae 
of the generalized plant of Fig. 12 and the formulae of the 
control object model that is part of the generalized plant, 
as described in the expressions (1) to (5), (11), and (32) 
to (40) . The control variable weight W yL for the L-th 
controlled variable y L is set in the control variable weight 
input unit 7 by the user of the design device. This control 
variable weight W yL is set for each controlled variable y. 
The control variable weight registration unit 8 outputs the 
control variable weight W yL input from the control variable 
weight input unit 7 to the weight matrix calculation unit 
9. The weight matrix calculation unit 9 calculates the weight 
matrix S, based on the control variable weight W yL , employing 
the expressions (29) and (30) , and outputs it to a controller 
calculation unit 6a. The controller calculation unit 6a 
registers the scaling matrix T and the weight matrix S in 
the formulae of the generalized plant stored in the memory 
unit 3a, and through y iteration, calculates the parameters 
of the controller K. At this time, the complementary 
sensitivity weight W t , the sensitivity weight W 3 and the 
scaling matrix M are preset in the generalized plant in the 



memory unit 3a . In this way, the controller K can be designed. 

As described above, in this third embodiment, each 
controlled variable y can be directly weighted by introducing 
the weight matrix S. Thereby, the controller with higher 
control performance and enhanced stability can be designed. 
Also, there is no need of providing the scaling matrix T with 
a role of weighting the controlled variable y by introducing 
the weight matrix S. 

Fourth embodiment 

In the third embodiment, the control variable weight 
adjusting means 20 (weight matrix S) is provided inside the 
closed loop system, but may be provided outside the closed 
loop system. Fig. 14 is a block diagram showing the 
configuration of the generalized plant in a fourth embodiment 
of the present invention. In this fourth embodiment, the 
control variable weight adjusting means 20 (weight matrix 
S) is provided at the former stage of the frequency 
sensitivity weight adjusting means 13. A method for 
determining the weight matrix S is exactly the same as the 
method for determining in the third embodiment as described 

in the expressions (29) and (30) . 

Next, in this fourth embodiment, the configuration of 

the generalized plant is modified as shown in Fig. 14, so 

that the expression (4) can be rewritten into the following 

expression. 

z , ' =S e , ... (4 1) 

Thereby, the expression (9) can be rewritten in the 



following expression. 
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Representing the expressions (8), (10) and (42) in the 
Doile's notation as in the expression (11), the parameters 
Ci and C 2 in the expression (11) can be represented in the 
following expression. 



Ci = [ SD £ Cp S Q a ] - C 2 = LD a C p C a J • • • (4 3) 

10 Also, the parameters Dn, D 12 , D 2 i and D 22 can be 

represented in the following expression. 

Dn = |-sD a SDa D pl Mj Di2 = js D „D p2 Tj 

D J1 = | r D a D a D pl Mj , Dz,=D a D p2 T • • - (4 4) 

The parameters A, Bi and B 2 are the same as shown in the 
15 expressions (12) and (13) . In this fourth embodiment, the 
actual controller mounted in the plants such as distillation 
tower is the controller K multiplied by the weight a _1 I and 
the scaling matrix T, as shown in Fig. 6. 

In this fourth embodiment, the configuration of the 
20 design device is almost the same as in the third embodiment. 
Thus, the operation of the design device of this fourth 
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embodiment will be described below with reference to Fig. 
11 . 

The operation of the control object model input unit 
1, the control object model registration unit 2, the frequency 
response calculation unit 4 and the scaling matrix 
calculation unit 5 is exactly the same as in the first or 
second embodiment. A memory unit 3a stores the formulae of 
the generalized plant of Fig. 14 and the formulae of the 
control object model that is part of the generalized plant, 
as described in the expressions (1) to (3), (5) to (8), (10) 
to (13), and (41) to (44). The operation of the control 
variable weight input unit 7, the control variable weight 
registration unit 8 and the weight matrix calculation unit 
9 is exactly the same as in the third embodiment. The 
5 controller calculation unit 6a registers the scaling matrix 
T and the weight matrix S in the formulae of the generalized 
plant stored in the memory unit 3a, and through y iteration, 
calculates the parameters of the controller K. In this way, 
the controller K can be designed. 
0 In this fourth embodiment, the control variable weight 

adjusting means 20 (weight matrix S) is provided at the former 
stage of the frequency sensitivity weight adjusting means 
13, but may be provided at the latter stage of the frequency 
sensitivity weight adjusting means 13. Since the magnitudes 
5 of the gains of the control object model 11 are made consistent 
by the scaling matrix T, the weight matrix S can be easily 
adjusted, and it is prerequisite that the method for 



determining the scaling matrix T as described in the first 
or second embodiment is employed in the third or fourth 
embodiment . 

Fifth embodiment 
5 Fig. 15 is a block diagram showing the configuration 

of a design device of a controller in a fifth embodiment of 
the present invention. The design device of Fig. 15 comprises 
a transient response parameter input unit 101 for inputting 
the transient response parameters indicating the transient 
10 response characteristics of a closed loop system composed 
of the control objects and the controller, a transient 
response parameter registration unit 102 for registering the 
transient response parameters within the device, a closed 
loop transfer function calculation unit 103 for calculating 
15 the transient response characteristics of the closed loop 
system based on the transient response parameters input from 
the transient response parameter registration unit 102, a 
frequency sensitivity weight calculation unit 104 for 
calculating the frequency sensitivity weight for determining 
20 the set value followup characteristic of the closed loop 
system, based on the transient response characteristics of 
the closed loop system, and a controller calculation unit 
105 for calculating the parameters of the controller by 
applying the frequency sensitivity weight to the preset 
25 generalized plants. The transient response parameter input 
unit 101, the transient response parameter registration unit 
102 and the closed loop transfer function calculation unit 
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103 constitute setting means for setting the transient 
response characteristics of the closed loop system. 

An algorithm for designing the controller in accordance 
with the H oo logic is designed based on the generalized plants 
as represented using the control objects. Therefore, the 
generalized plants will be first described. Fig. 16 is a block 
diagram showing the configuration of a model of an actual 
control object as represented in the formulae. A numerical 
model of the control object as shown in Fig. 16 is composed 
of a model 111 for the manipulated variable u and a model 
112 for the disturbance w. Reference sign Pu denotes a 
transfer function of the model 111 and reference sign Pw 
denotes a transfer function of the model 112 . The models 111, 
112 are obtained as a result of the model identification using 
the data obtained by a step response test for the actual 
control object. The controlled variable y that is the output 
of control object is a sum of the outputs from the models 
111 and 112. 

Fig. 17 shows the configuration of the conventional 
generalized plants containing such numerical model of 
control object. The generalized plants involve providing a 
frequency weight called a sensitivity weight W s used to 
determine the set value followup characteristic and a 
frequency weight called a complementary sensitivity weight 
W t used to determine the robust stability, and introducing 
the set value r, and the outputs zl, z2 in addition to the 
input (manipulated variable) u of control object, input 



(disturbance) w, and output (controlled variable) y, to 
realize both the set value followup characteristic and the 
robust stability as shown in Fig. 17. Deviation e (=y-r) is 
an observation quantity or the input into the controller (not 
5 shown) . Reference numeral 113 denotes a block representing 
the, sensitivity weight W s , and reference sign Z 1 denotes an 
output for effecting the evaluation of the set value followup 
characteristic. Also, reference numeral 114 denotes a block 
representing the complementary sensitivity weight W t , and 
10 reference sign Z 2 denotes an output for effecting the 
evaluation of the robust stability. 

Conventionally, in the generalized plants as shown in 
Fig. 17, the complementary sensitivity weight W t is 
determined by estimating the uncertainty of the models on 
15 the basis of the numerical models of control objects, the 
sensitivity weight W s is determined by directly specifying 
the frequency characteristic in view of the f ollowupability 
to the set value r, whereby the parameters of the controller 
are determined through the y iteration. However, since 
20 employing the generalized plants of Fig. 17, the controller 
is designed on the basis of a larger gain due to differences 
in the gains for the plant outputs of the manipulated 
variables, the obtained controller is likely to be very 
conservative, or excessively stable. Also, since the set 
25 value followup characteristic and the disturbance response 
characteristic are usually reciprocal, it is preferred to 
design with the weight according to the purpose rather than 
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with the same weight. Further, in the generalized plants of 
Fig. 17 containing no integral element, some steady-state 
deviation arises. If the sensitivity weight W s is provided 
with the integral characteristic, the controller can have 
the integral characteristic, but the generalized plants 
become unstable, not leading to the normal H oo problem. 

Therefore, in this fifth embodiment, the generalized 
plant as shown in Fig. 18 is considered. In Fig. 18, reference 
sign M denotes a scaling matrix for adjusting the influence 
of disturbance w on the controlled variable y, reference sign 
T denotes a scaling matrix for making the magnitudes of errors 
of the control objects consistent, and reference sign a' 1 ! 
denotes a weight for providing the controller with the 
integral characteristic to eliminate the steady-state 
deviation. Herein, a(s) is defined as a(s)= s/ (s+a) . Where 
s is a Laplace operator, and a (>0) is any real number. 
Reference numeral 115 denotes a block representing the 
scaling matrix M, reference numeral 116 denotes a block 
representing the scaling matrix T , and reference numeral 117 
0 denotes a block representing the weight a _1 I. Deviation e 2 
is deviation e multiplied by the weight a -1 I, and input into 
the controller. Fig. 19 shows the configuration of a robust 
control system in which a controller K is added to the above 
generalized plant. In Fig. 19, reference numeral 118 denotes 
5 a block representing the controller K. 

The design device of the controller in this fifth 
embodiment is aimed at determining the parameters of the 



controller K such that the controlled variable y that is the 
output of control object follows the set value r, the 
influence of disturbance w is removed, and the control object 
is stabilized even if it is fluctuated or there is an error 
in the model of control object. The H oo control problem can 
be regarded as the problem of reducing the H oo norm (gain) 
of the transfer function from (r, w) to [z lf z 2 ) . That is, 
the set value followup characteristic, the robust stability 
and the disturbance suppression may be considered in the 
following way. 

(D) Set value followup characteristic: If the H oo norm 
(gain) of the transfer function from the set value r to 
deviation e (more correctly, transfer function with the set 
value r multiplied by the frequency weight a _1 W s from r to 
zi) is reduced, the deviation e can be decreased, so that the 
set value followup characteristic can be made better . Herein, 
a _1 W s is the frequency weight for restricting the followup 
band (e.g., followup only in the low band) . 

(E) Robust stability: There is an error between the 
actual control object and its model due to characteristic 
variations of the control object or the error at the time 
of modeling. The maximum value of the error from the 
identified model is estimated as A(s), and the controller 
K is designed such that the H oo norm from the set value r 
to z 2 is smaller than or equal to 1, employing the 
complementary sensitivity weight W t (s) such that 
I A (jo) |< I W t (jo) I for this error, whereby the robust 



stabilization can be achieved. 

(F) Disturbance suppression: If the H oo norm (gain) of 
the transfer function from disturbance w to deviation e (more 
correctly, transfer function with the disturbance w 
multiplied by the frequency weight cc _1 W s from w to z 1 ) is 
reduced, the deviation e can be decreased even if the 
disturbance w enters, whereby the disturbance suppression 
can be ameliorated. 

Next, it is supposed that the state space representation 
of the generalized plant as shown in Fig. 18 is given by: 

X p =A p x p + B pl Mw+B p2 Tu •■• (101> 

y = C p x p +D pl Mw + D p2 Tu (102) 

In the above expressions (101) and (102) , x p is the 
quantity of state, and A p , B pl , B p2 , C p , D pl and D p2 are the 
parameters of the numerical models 111, 112 of control objects. 
From the expression (102), the deviation e can be obtained 
in the following expression. 

e = y-r = C p x p +D pl Mw + D p2 Tu-r •** (103) 

With the configuration of the generalized plant as shown 
in Fig. 18, the outputs z a ' and z 2 ' can be defined in the 
following expressions. 

z,' =e, * ' " (104) 

Z2 . =u - - • (1 0 5) 

The frequency weight for providing the controller K with 

the integral characteristic can be defined in the following 



expression, using the expression (103) . 

x a = A a x a + B a e 

= A a x a + B a C p x p + B a D pl M w + B a D p2 T u - B a r •••(106) 



e 2 = C a x a + D a e 

= C a x a + D a C p x p + D a D pl M w 4- D a D p2 T u - D a r •••(107) 

In the expressions (106) and (107), x a is the quantity 
of state of a" 1 1, and A«, B a , C a and D a are the parameters of 
a -1 1. Arranging the above expressions and representing them 
in a state space, the following three expressions can be 
obtained. 



[ B <A> A «][ X a] + [- B ' 



:;H D o Cp c oi[ x x :H~o° D «r Da : p2 1 

=2 = |P«Cp C a ^]+[-D a D a D pl M D a D p2 T|w 



(10 9) 



(1 1 0) 



Representing the expressions (108), (109) and (110) in 
Doile's notation, the following expression can be obtained. 



G(s)= 



A B x 
Ci D u 
Co D a 



B 2 ' 
D 12 



(1 1 1) 



Parameter A can be represented as follows: 





■ 0 



B p iM ■ 
B a D pl M 




(1 1 3) 



Also, parameters d and C 2 can be represented as follows : 



And the parameters Du, D 12 , D 2 i, D 22 can be represented 
as follows: 



The sensitivity weight W s and the complementary 
sensitivity weight W t are designed, and multiplied by the 
output parts of the expression (111) . Through y iteration, 
the controller K is obtained in the state space representation. 
Herein, the output parts of the expression (111) signify the 
parts corresponding to the outputs zi' , z 2 ' in Fig. 18. Hence, 
an output equation of the parameters Ci, D 1X , Di 2 of the 
expression (111) may be multiplied by a diagonal matrix Q 
as represented in the following expression having the 



(1 1 4) 




(1 1 5) 



diagonal elements of sensitivity weight W s and complementary 
sensitivity weight W t from the left side. Thus, the parameters 
of the controller K can be calculated. 

Q= [ W * °1 •••(116) 

The controller K is a solution of the H oo control problem 
with the generalized plants, and the actual controller 
mounted on the plants such as a distillation tower is the 
controller K multiplied by weight a 1 ! and scaling matrix T, 
as shown in Fig. 20. 

As previously described, the H oo logic is a design method 
in the frequency domain. Therefore, the design is easy in 
the control of the mechanical system, but is hard in the 
control system untreatable in the frequency domain such as 
the process control, and it is difficult to suitably select 
the complementary sensitivity weight W t and the sensitivity 
weight W s . In the following, a method for determining the 
complementary sensitivity weight W t in this embodiment will 
be described below. The control object is varied in the 
characteristics, depending on the driving conditions. 
Normally, the control design is made based on a certain model, 
but in the robust control design, the variation of control 
object and the magnitude of error of modeling are contained 
beforehand in the control design, and the controller is 
designed to be stable without much deterioration in the 
control performance even if there is any variation or error. 



Fig. 21 shows an additive error for the model 111 of control 
object. In Fig. 21, reference numeral 119 denotes a block 
representing the additive error A. In the robust control 
design, a variation in the characteristic of the control 
5 object owing to the driving conditions and a model error due 
to lower dimension of the model 111 are represented as the 
additive error A as shown in Fig. 21. If the characteristics 
of the control object are deviated from the model 111 due 
to this additive error A, the controller is designed such 
10 that the controller output may be stable. To this end, the 
complementary sensitivity weight W t may be determined to 
cover the additive error A. The general expression of this 
complementary sensitivity weight W t is shown in the following 
expression. Since the change of the model 112 is not related 
15 with the stability of the system, it is supposed that the 
model 111 alone is varied in designing the controller. 
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In this fifth embodiment, for the additive error A with 
20 scaled magnitude, employing the scaling matrix T, the element 
of the complementary sensitivity weight W t is the maximum 
value Gmax of the gains of error A multiplied by a safety 
factor 5 (5 is equal to 1, for example) . Namely, the elements 
(weights) W t i, W t2 , W t3 , W tN are defined as follows: 

- 49 - 



W ll = W l2 = W, 3 = W t N= (1 + 6) Gmax • • • (118) 

This embodiment is involved with the multivariable 
control system, and assuming that the number of manipulated 
variables u is N (N is a positive integer) , the complementary 
sensitivity weight W t is NxN matrix. Where W tN is the weight 
for the N-th manipulated variable u N . 

Next, the methods for determining the scaling matrices 
T, M will be described below. A general expression of the 
scaling matrix T is shown in the following' expression. 



T x 0 0 

0 T 2 0 

0 0 T 3 

0 0 0 



(119) 



Assuming that the number of manipulated variables u is 
N, the scaling matrix T is NxN matrix. Element T N of the scaling 
matrix T is the weight for the N-th manipulated variable u N . 
Each element T N is determined such that the magnitudes of the 
components of the additive error A are as equal as possible. 

Next, a general expression of the scaling matrix M is 
shown in the following expression. 



M = 



M x 0 0 

0 M 2 0 

0 0 M 3 

0 0 0 



(12 0) 



Assuming that the number of disturbance w is J (J is 



a positive integer), the scaling matrix M is JxJ matrix. 
Element M a of the scaling matrix M is the weight for the J-th 
disturbance Wj. Each element Mj is an adjustment parameter 
for determining the disturbance suppression performance by 
5 adjusting the influence of each disturbance Wj on the 
controlled variable y. 

Next, a method for determining the sensitivity weight 
W s in this fifth embodiment will be described below. First 
of all, a general expression of the sensitivity weight W s is 
10 shown in the following expression. 
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Assuming that the number of controlled variables y is 
L (L is a positive integer) , the sensitivity weight W s is LxL 

15 matrix. Element W sL of the sensitivity weight W 3 is the weight 
for the L-th controlled variable y L . In order to determine 
the sensitivity weight W s , a closed loop system of Fig. 22 
in which the robust control system of Fig. 19 is simplified 
is considered. In Fig. 22, reference numeral 111a denotes 

20 a block representing a numerical model P of control object, 
and reference numeral 113a denotes a block representing the 
frequency weight W s ' . 

Supposing that the sensitivity function indicating the 
control performance mainly regarding the quick-response 
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property such as the set value followupability or disturbance 
suppression is S(s), the smaller gain |S(j©)| of the 
sensitivity function S(s) is preferable because the model 
variations have less effect on the set value response. If 
the control specification at each frequency is given by 
S 3P ec(©), the following condition concerning the sensitivity 
function S(s) can be obtained. 

|S(j©)| < S sp ec(a>); Vg> (122) 

Where V© means that the expression (122) holds for all 
frequencies ©. 

On one hand, the sensitivity function S(s) corresponds 
to the transfer function G er (s) of the closed loop system from 
the set value r to deviation e, as shown in Fig. 22. Supposing 
that the transfer function of control object is P(s) and the 
transfer function of controller is K(s), the following 
expression is obtained. 



s W- G ««-rr^) •■• d2 3) 

Also, the transfer function G yr (s) of the closed loop 
system of Fig. 22 from the set value r to the controlled 
variable y is obtained in the following expression. 



G ( s )= " P ^A • • • (12 4) 

y A) l-P(s)K(s) 

From the expressions (123) and (124), the sensitivity 

function S (s) can be obtained in the following expression. 



- 52 - 



S(s)=G yr (s)-l • ' ' (125) 

In this fifth embodiment, the transient response 
characteristic of the closed loop system as shown in Fig. 
22, or the transfer function G yr (s) of the closed loop system 
from the- set value r to the controlled variable y is 
approximated with a first-order lag characteristic as in the 
following expression, giving G yr - spec (s) as the specification 
of the transfer function G yr (s). 



G ( s }= * • • • (12 6) 

In the expression (126) , T sL is a time constant regarding 
the L-th controlled variable y L . Employing G yr - spec (s) in the 
expression (126), instead of G yr (s) in the expression (125), 
the following control specification S spe c(s) is obtained. 



s ( s )= - T *L S ... (127) 

* Sp ecKV T sL S + 1 

In this fifth embodiment, the frequency weight W sL ' (s) 
regarding the L-th controlled variable y L is set as follows: 

'-^-sdsj"-^ ••• (128> 

The frequency weight W sl ' (s) is W sL (s) multiplied by 
a _1 (s), and defined as follows: 
W. L '«=a-^)W BL W ' • ' (12 9) 
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If the frequency weight W sL ' (s) is set as in the 
expression (128), the expression (122) can be transformed 
into the following expression. 



Further, the expression (130) can be rewritten into the 
following expression, using the H oo norm. 

|w sL (s)S(s)J 0; <l • • • (13 1) 

The expression (131) indicates that the H oo norm of the 
transfer function of the closed loop system of Fig. 22 from 
the set value r to the deviation e (more correctly, the 
transfer function with the set value r multiplied by the 
frequency weight a _1 (s)W sL (s) from r to z x ) is less than 1. 
This expression (131) is a design index of the controller 
K in view of the set value followup characteristic. 
Accordingly, by setting the frequency weight W sL ' (s) in 
accordance with the expression (128), the expression (131) 
is satisfied, whereby the controller K can be designed in 
view of the set value followup characteristic. If the 
expression (128) is transformed, the following expression 
is obtained. 

^'^•W^T (132) 

The first term on the right side of the expression (132) 
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isa _1 (s) . Therefore, element W sl (s) of the sensitivity weight 
W s can be calculated as follows: 



W sL (s) = 



T s r s + 1 



(13 3) 



T sL (s + a) 

Substituting the expression (133), the expression (121) 
can be represented as follows: 



T sl (s+a) 
0 

0 



0 



_ T s2 s + 1 
T s2 (s+a) 



0 



. T s3 s + 1 
T s3 (s+a) 



0 



0 
0 
0 

T s lS + 1 



T sL (s + a) 
• • • (13 4) 

In this way, the sensitivity weight W s can be determined. 
Referring now to Fig. 15, the above operation will be 
described below. The transient response parameter or the time 
constant T sL is set in the transient response parameter input 
unit 101 by the user of the design device. The time constant 
T sL is set for each controlled variable y. The transient 
response parameter registration unit 102 outputs the time 
constant T sL input from the transient response parameter, input 
unit 101 directly to the closed loop transfer function 
calculation unit 103. The closed loop transfer function 
calculation unit 103 calculates the transfer function G yr (s) 
of the closed loop system as shown in Fig. 22 from the set 
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value r to the controlled variable y by substituting the input 
time constant T sL into the expression (126) , and outputs it 
to the frequency sensitivity weight calculation unit 104. 
Subsequently, the frequency sensitivity weight calculation 
unit 104 calculates the sensitivity weight W s , based on the 
transfer function G yr (s), employing the expressions (125), 
(127), (128), and (132) to (134), and outputs it to the 
controller calculation unit 105. The memory unit 106 stores 
the formulae of generalized plants as described in the 
expressions (101) to (115) and shown in Fig. 18. The 
controller calculation unit 105 registers the sensitivity 
weight W s in the formulae of the generalized plants stored 
in the memory unit 106, and through y iteration, calculates 
the parameters of the controller K. At this time, the 
complementary sensitivity weight W t , and the scaling matrices 
T, M are preset in the generalized plants in the memory unit 
106. In this way, the controller K can be designed. 

As previously described, in the process control field, 
it is difficult to give the control specification in the 
frequency domain, and to determine the frequency sensitivity 
weight W s . On the contrary, in this embodiment, the transient 
response characteristic (transfer function G yr (s)) of the 
closed loop system is set, and the frequency weight W s can 
be calculated, based on this transient response 
characteristic. Thereby, even in the process control field, 
where the frequency response characteristic is difficult to 
give as the control specification, the controller can be 



designed in accordance with the H oo logic. Consequently, the 
multivariable control system can be easily designed in view 
of the variations of the control objects and the uncertainty 
of the numerical models. Also, the controller can be realized 
utilizing the features of the H <x> logic with smaller 
computational load during execution of the control, and with 
which the small-scale control system can be implemented. 
By approximating the transient response characteristic of 
the closed loop system with a first-order lag characteristic, 
the parameters for the control design are intuitively 
understandable to the designer. Therefore, the design device 
that is easily understood and used by the designer can be 
realized. Since the parameters are intuitively 
understandable to the designer, the design device that is 
capable of changing the design can be realized in the case 
where the controller is redesigned after designed once. 
Sixth embodiment 

In the fifth embodiment, the transient response 
characteristic of the closed loop system is approximated with 
a first-order lag, and the sensitivity weight W s is determined 
by specifying the time constant of the closed loop system 
for the controlled variables . In this case, albeit the control 
design based on the transient response characteristic, only 
one design parameter is specified for each controlled 
variable, and the obtained control system does not 
necessarily operate with the first-order lag as specified. 
Considering the practicality of the multivariable control, 
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it is not favorable to increase the number of design 
parameters, but if the degree of freedom for the design is 
increased by making the meanings of design parameters clear, 
it may not be necessarily a false policy. Thus, in this sixth 
embodiment, the controller is designed by approximating the 
characteristics of the closed loop system with a second- 
order system as commonly used, and specifying the transient 
response characteristic parameters. 

Firstly, the transient response of the closed loop 
system of Fig. 22 for each controlled variable is approximated 
with the second-order system in the following expression. 

y- W n 2 ... (135) 

r s 2 + 2 C w n s + co n 2 

In the expression (135), C (C>°) is a damping coefficient , 
and con (C>°) is a natural frequency. The solution of the 
characteristic equation of the system as indicated in the 
expression (135) is obtained in the following expression. 



s=-Cco n ±o>Jc 2 -l • • • (136) 

Giving a step input (r (s) = 1/s) , the controlled variable 
y(s) that is an output of the closed loop system can be 
obtained in the following expression. 



y(s) (s- Sl Xs-s 2 )s 

= 1 + C i ° 2 • • • (137) 

s s + Co) n -a) n ^ 2 -l s + £co n +co n V£ 2 -l 
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The parameters Si, s 2 , Ci, C 2 in the expression (137) are 
given in the following expression. 



Sl = -Co) n +co nA /£ 2 -l . s 2 =-Cco n -co n> f? 



C 7 = ^ • • - (13 8) 



Consider the most important damped oscillation (0 < C 
< 1) in the second-order system characteristic. Then, the 
parameters Ci, C 2 are given in the following expression. 



Cl = _i / I , C 2 =-A J • • • (139) 

Thereby, the time response of the closed loop system 
as shown in Fig. 22 can be represented in the following 
expression. 
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y ( t )=l+_l / -exp -C^ + i coJl-C 2 U 

2V1-C 2 V ' 



if-C + iVl-C 2 ) , 

+ _i i-expf - £ con - i <o n V 1 ~ £ * 

2V1 -c 2 V ' 

= 1 + e ~ c<0nt f cosco n Vl - S 2 t + isitico n ^l - C 2 t) 

2 A /T :r c 2 v / 
+ iJi-cM , , , — \ 

+ i == L. e-^^fcos^Vl^^t- isinca„Vl-C 2 t ) 

2V1- £ 2 ^ ' 

= l- e - Ct0nt |cos6) ll> /l-C 2 t+ . C sineo n> /l-C 2 t 

I Vi-c 2 J 

_ 1 _ e~ c<0nt ^/ ^ ^ _ y 2 f4 .^ > ^cos^C - ■ • (140) 

One embodiment of the time response of the closed loop 
system as represented in the expression (140) is shown in 
Fig. 23. Fig. 23 shows how the controlled variable y behaves 
when a 100% step set value r is given at time 0. The parameters 
representing the typical transient response characteristic 
of the second-order system corresponding to this time 
response include a rise time t r for which the controlled 
variable y reaches the same value as the set value r (herein 
100%) , the overshoot OS that is an extreme value of the 
transient deviation firstly taken after the controlled 
variable y exceeds the set value r, the overshoot time t p for 
which the controlled variable y reaches the overshoot OS, 
the settling time t q for which the controlled variable y is 
contained within 5% of a range of the set value r, and the 
damping ratio DR that is the ratio of a x to a 2 as shown in 
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Fig. 23. The rise time t r can be represented in the following 
expression, employing the damping coefficient C an d the 
natural frequency co n . 

t r = J—^TC-cos-^) ' ' • (14 1) 

a>Jl-Z 2 

Similarly, the overshoot time t p is indicated below. 

••• (142) 



P 

The settling time t q can be obtained in the following 
expression. 

t q =7^- • ' • (143) 

Also, the overshoot OS can be obtained in the following 
expression, employing the damping coefficient 



And the damping ratio DR can be obtained in the following 
expression . 



DR = a 2 /a 1 =exp^-2 TcC/Vl-^ 2 ) *** < 145) 

From the expression (135) , the transfer function G yr (s) 
of the closed loop system of Fig. 22 from the set value r 
to the controlled variable y can be obtained in the following 
expression . 

G yr (s)=- ^ ... (146) 
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If the frequency weight W sL ' (s) is set in the following 
expression, as in the fifth embodiment, the expression (131) 
is satisfied, and the controller K can be designed in view 
of the set value followup characteristic. 

> M _ 1 _ 1 = 1 

WsL(S) ~V^)~G yr (s)-l co nL 2 



S 2 + 2C L o) nL s+co nL 2 

• • • (14 7) 

In the expression (147), Cl is the damping coefficient 
for the L-th controlled variable y L , and ra nL is the natural 
frequency for the controlled variable y L . The expression 
(147) can be transformed into the following expression. 

w . sii^s^ 
- s 2 -2 Cl^itLS 
s + a -(s 2 + 2£ T Q) nL s+co nL 2 ) _ m < (148) 

s (s + 2 C L o) nL Xs + a) 

The first term on the right side of the expression (148) 
is a _1 (s) . Accordingly, the element W sL (s) of the sensitivity 
weight W s can be calculated in the following expression. 

W M s 2 + 2 ^n L s+co nL 2 # a m (149) 

Wsl(s) " (s + 2£ L co nL Xs + a) U49) 

Substituting the expression (149), the expression (121) 

can be represented as follows: 
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' s 2 +2C 1 a> nl s-hfa) nl 2 
(s + ^mXs + a) 



S 2 +2C 2 C0 n2 S + 0J n2 2 

(s+2C 2 co n2 Xs+a) 



0 



0 



0 



0 



0 



0 



S 2 +2C L C0 nL S + 6> nL 2 



(s + 2S L w„LX s + a ) 



(15 0) 



To calculate the sensitivity weight W s employing the 
expression (150), it is required to obtain the damping 
coefficient Cl and the natural frequency © n L- To obtain the 
5 damping coefficient Cl and the natural frequency ©nL/ two items 
are selected from among the rise time t r , the overshoot time 
t p , the settling time t q , the overshoot OS and the damping 
ratio DR, and the two selected parameter values are set. In 
this sixth embodiment, the method for calculating the damping 
10 coefficient Cl and the natural frequency aw as one embodiment, 
will be described with an instance of employing the rise time 
t r and the overshoot OS. With the overshoot OS as indicated 
in the expression (144), the damping coefficient Cl can be 
obtained in the following expression. 



In the expression (151), OS L is the overshoot for the 
L-th controlled variable y L . Also, from the rise time t r as 
indicated in the expression (141) , the natural frequency © nL 
can be obtained in the following expression. 



15 




(15 1) 
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In the expression (152), t rl is the rise time for the 
L-th controlled variable y L . In this way, if the value of 
overshoot 0S L is set, the damping coefficient Cl can be 
5 calculated from the expression (151), and if the value of 
the rise time t rL is set, the natural frequency (D nL can be 
calculated from the expression (152) . In the above way, the 
sensitivity weight W s can be determined. 

In this sixth embodiment, the configuration of the 
10 design device is most the same as in the fifth embodiment. 
Thus, the operation of the design device of this sixth 
embodiment will be described below with reference to Fig. 
15. The transient response parameters, i.e., the overshoot 
OS L and the rise time t rL , are set into the transient response 
15 parameter input unit 101 by the user of the design device. 
This overshoot OS L and the rise time t rL are set for each 
controlled variable y. The transient response parameter 
registration unit 102 outputs the overshoot OS L and the rise 
time t rL input from transient response parameter input unit 
20 101 directly to the closed loop transfer function calculation 
unit 103. The closed loop transfer function calculation unit 
103 calculates the damping coefficient Cl and the natural 
frequency co nL , based on the overshoot OS L and the rise time 
t rL / employing the expressions (151) and (152) . Substituting 
25 the damping coefficient Cl and the natural frequency © nL into 
the expression (146), the transfer function G yr (s) is 
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calculated, and output to the frequency sensitivity weight 
calculation unit 104. Then, the frequency sensitivity weight 
calculation unit 104 calculates the sensitivity weight W a , 
based on the transfer function G yr (s), employing the 
expressions (147) to (150), and outputs it to the controller 
calculation unit 105. The operation of the controller 
calculation unit 105 and the memory unit 106 is exactly the 
same as in the fifth embodiment. In this way, the controller 
K can be designed. 

As described above, in this sixth embodiment, the degree 
of freedom in the design can be increased by approximating 
the transient response characteristic of the closed loop 
system with the second-order system characteristic, and the 
controller provided by the design device can have wider 
applicability. The first to sixth embodiments are involved 
with the design device which designs the multivariable 
controller. Also, the design device of the first to sixth 
embodiments can be implemented on the computer. That is, the 
computer is equipped with an operation unit, a storage device 
and an input/output device, and operates as the design device 
in accordance with the program. 
Industrial Applicability 

As described above, the present invention is suitable 
for designing the multivariable controller. 
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CLAIMS 

1. A design device for designing a controller in 
accordance with the H infinity (H co) control logic, employing 
the generalized plants having the control object models for 
the manipulated variables, characterized by comprising: 
storage means for storing said generalized plants; 
parameter calculating means for calculating the 
parameters of components of said generalized plants in 
accordance with the response characteristic of said control 
object model or the response characteristic of a closed loop 
system consisting of said control object model and said 
controller; and 

controller calculation means for calculating said 
controller by applying said parameters to said generalized 
plants stored in said storage means. 

2. The design device of controller according to claim 
1, characterized in that said generalized plants have said 
control object model, and manipulated variable weight 
adjusting means for adjusting the input of manipulated 
variable into said control object model, which is provided 
in the former stage of said control object model, said 
parameter calculating means comprises frequency response 
calculation means for calculating the frequency responses 
of said control object models, and scaling matrix calculation 
means for calculating a scaling matrix T for determining the 
weighting of said manipulated variables with said 
manipulated variable weight adjusting means in accordance 



with the frequency responses of said control object models 
so that the respective gains of said control object models 
are consistent, and said controller calculation means 
calculates the controller by applying said scaling matrix 
T to the manipulated variable weight adjusting means of said 
generalized plants stored in said storage means. 

3. The design device of controller according to claim 
2, characterized in that said scaling matrix calculation 
means calculates said scaling matrix T as follows: 



T = 



r x o 

0 To 



0 0 T 3 



_1 mtt\|Gyiui|, 



0 
0 
0 



yiu2j| ( 



-l G : 



yluN||, 



J yiuN 



^ x max (||Gy 2m | oo , | G y2u2|| 0O > - » l G y2u N l J 
L || G y2uN|| 00 

+ i max (iGymiL » l^y^L > - > ISlunL ) 

where the number of manipulated variables u is N (N is a 
positive integer) , the number of controlled variables y is 
L (L is a positive integer) , and the H qo norm of the transfer 
function of said control object model from the N-th 
manipulated variable u N to the L-th controlled variable y L 

iS ||G yLu N||oo. 

4. The design device of controller according to claim 
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1, characterized in that said generalized plants have a first 
control object model for the manipulated variables, a second 
control object model for the disturbance, and manipulated 
variable weight adjusting means for adjusting the input of 
manipulated variable into the first control object mode, 
which is provided in the former stage of said first control 
object model, said parameter calculating means comprises 
frequency response calculation means for calculating the 
frequency responses of said first control object model and 
said second control object model, and scaling matrix 
calculation means for calculation a scaling matrix T for 
determining the weighting of the manipulated variables with 
said manipulated variable weight adjusting means in 
accordance with the frequency responses of said first and 
second control object models so that the respective gains 
of said first control object model are consistent with the 
maximum values of the gains of said second control object 
model, and said controller calculation means calculates the 
parameters of said controller by applying said scaling matrix 
T to the manipulated variable weight adjusting means of said 
generalized plants stored in said storage means. 

5. The design device of controller according to claim 
4, characterized in that said scaling matrix calculation 
means calculates said scaling matrix T as follows: 
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T = 



■T t 0 0 
0 T 2 0 
0 0 T 3 



0 
0 
0 



0 



0 



0 




1 max 



(|Gyiwi|.-l G yiw2L'--l G yw4) 





max 

L* 



I G y2uN|L 

™a* (| , 1 , - , | Gy^ ) 




where the number of manipulated variables u is N (N is a 
positive integer) , the number of disturbance w is J (J is 
a positive integer) , the number of controlled variables y 
5 is L (L is a positive integer) , the H oo norm of the transfer 
function of said first control object model from the N-th 
manipulated variable u N to the L-th controlled variable y L 
is ||G yL un||oo, and the H oo norm of the transfer function of said 
second control object model from the J-th disturbance w a to 
10 the L-th controlled variable y L is ||G yLw j|U- 

6. The design device of controller according to claim 
2, characterized in that the generalized plants stored in 
said storage means have control variable weight adjusting 
means for adjusting the controlled variables inside a closed 
15 loop system consisting of said manipulated variable weight 
adjusting means, the control object model and the controller, 
and said design device has setting means for setting a weight 
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matrix S for determining the weighting of the controlled 
variables with said control variable weight adjusting means . 

7. The design device of controller according to claim 
4, characterized in that the generalized plants stored in 
said storage means have control variable weight adjusting 
means for adjusting the controlled variables inside a closed 
loop system consisting of said manipulated variable weight 
adjusting means, the first control object model and the 
controller, and said design device has setting means for 
setting a weight matrix S for determining the weighting of 
the controlled variables with said control variable weight 
adjusting means. 

8. The design device of controller according to claim 
2, characterized in that the generalized plants stored in 
said storage means have control variable weight adjusting 
means for adjusting the controlled variables in the former 
or latter stage of frequency sensitivity weight adjusting 
means for determining the set value followup characteristic 
of a closed loop system consisting of said manipulated 
variable weight adjusting means, the control object model 
and the controller, and said design device has setting means 
for setting a weight matrix S for determining the weighting 
of the controlled variables with said control variable weight 
adjusting means. 

9. The design device of controller according to claim 
4, characterized in that the generalized plants stored in 
said storage means have control variable weight adjusting 



means for adjusting the controlled variables in the former 
or latter stage of frequency sensitivity weight adjusting 
means for determining the set value followup characteristic 
of a closed loop system consisting of said manipulated 
5 variable weight adjusting means, the first control object 
model and the controller, and said design device has setting 
means for setting a weight matrix S for determining the 
weighting of the controlled variables with said control 
variable weight adjusting means. 
10 10. The design device of controller according to claim 

1, characterized in that said parameter calculating means 
has setting means for setting the transient response 
characteristic of said closed loop system and frequency 
sensitivity weight calculation means for calculating the 
15 frequency sensitivity weight for determining the set value 
followup characteristic of said closed loop system in 
accordance with the transient response characteristic of 
said closed loop system, and said controller calculating 
means calculates the parameters of said controller by 
20 applying said frequency sensitivity weight to said 
generalized plants stored in said storage means. 

11. The design device of controller according to claim 
10, characterized in that said frequency sensitivity weight 
calculation unit calculates said frequency sensitivity 
25 weight in accordance with the transient response 

characteristic of said closed loop system, and a design index 
that the H oo norm of the transfer function of the closed loop 
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system from the set value to the deviation multiplied by said 
frequency sensitivity weight is less than 1. 

12. The design device of controller according to claim 
10, characterized in that said setting means approximates 
the transient response characteristic of said closed loop 
system with a first-order lag characteristic. 

13. The design device of controller according to claim 
10, characterized in that said setting means approximates 
the transient response characteristic of said closed loop 
system with a second-order system characteristic. 
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